U. S. Department of Commerce National Bureau of Standards 

RESEARCH PAPER RP1493 

Part of Journal of Research of the Rational Bureau of Standards, Volume 29, 

August 1942 

AN EQUATION FOR THE ISOTHERMS OF PURE SUB- 
STANCES AT THEIR CRITICAL TEMPERATURES 

By Cyril H. Meyers 



abstract 

An'empiricaFequation^is presented for data~along the isotherms of pure sub- 
stances at their critical temperatures and densities up to 1.1 times the critical 
density. The equation has five constants in addition to the universal gas con- 
stant. Of these five, one is limited to a single value, if a simple solution to the 
equation is to be obtained; three are determined by the critical conditions; and 
the remaining constant is determined empirically from data other than those for 
the critical state. Three of the five constants are independent of the substance, 
one is an explicit function of RT c lp c V e , whereas the reciprocal of the fifth is found 
to be practically a linear function of RT c /p c V e for values of that ratio up to 3.8 or 
3.9. This relation correlates the isotherms for various substances, so that data 
on two well chosen substances, such as hydrogen and carbon dioxide, suffice to 
determine the fifth constant and consequently the critical isotherm for other sub- 
stances, for which RT c /p c V e is not greater than 3.8. For substances such as 
ammonia (for which RT c /p c V c is 4.08), the fifth constant may not conform to the 
linear relation. 

The equation represents the data for 12 of 13 substances within the experi- 
mental error, the exceptional substance being water (for which RT c /p c V e is 4.34). 

At low densities the equation reduces to a simpler form which expresses pV/RT 
as a quadratic function of density. «■';■ ., 

The fugacity along the critical isotherm can be calculated frorcfthe equation. 
The fugacity at the critical state is approximately two-thirds the critical pressure 
for all substances. The fungacity of CO2 is calculated at 1 atmosphere and at the 
critical pressure. 
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I. INTRODUCTION 

The isotherm at the critical temperature forms a convenient base 
from which to proceed in tabulating the properties of vapors. Conse- 
quently, a mathematical representation of the critical isotherm is 
desirable. A mathematical representation of the critical isotherm is 
especially difficult, since at the critical temperature not only must the 
equation represent the experimental data, but it must satisfy the 
conditions that the first and second derivatives of pressure with respect 
to volume (or density) must equal zero at the critical state. 

The notation used in this paper is as follows: R, the universal gas 
constant; T, absolute temperature; p, absolute pressure; V, molal 
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Figure 1. — Comparison of reduced pressures at corresponding densities along the 
critical isotherms of four substances. 

volume; c, subscript representing critical state; X, VJV or reduced 
density;/, fugacity; m, number of a term in a series; n, empirical 
constant; and K, empirical constant. 

II. FAILURE OF THE LAW OF CORRESPONDING STATES 

The law of corresponding states has been used occasionally for 
predicting the properties of pure substances from the properties of 
other substances, but this procedure is subject to rather large errors. 
According to the law of corresponding states, the data for the critical 
isotherms of all substances would lie on a single curve on a graph with 
reduced pressure as ordinate and reduced density (or volume) as 
abscissa. Such a graph, which is illustrated in figure 1 for H 2 , C0 2 , 
NH 3 , and H 2 0, indicates that the assumption of the law of corres- 
ponding states may lead to errors as large as 10 percent of the pressure. 
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For densities less than the critical density, the position of the curves 
appears to be related to the quantity RT c /p c V c . The greater the value 
of RTJpcVc, the greater is the reduced pressure corresponding to any 
given reduced density, X. This has been found to be true also for 
other substances studied but is omitted from figure 1 to avoid con- 
fusion. 

The arrangement for densities greater than the critical is in general 
reversed from that at smaller densities and appears not to be entirely 
a function of RT c /p c V c , although this irregularity may be due in part 
to erroneous critical densities. The data for water are so conspicuously 
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V c /V = X 

-Comparison of data for the critical isotherms of four substances. 



far below the other curves in figure 1 that the irregularity appears to 
be greater than any experimental error possible in the critical density. 

III. EQUATION FOR THE ISOTHERM AT THE CRITICAL 

TEMPERATURE 

The differences between the critical isotherms of different substances 
at densities less than the critical are better illustrated by plotting the 
quantity 

V(RT-pV) X 
V c (RT-p c V c y 2 

as ordinate, and the reduced density, X, as abscissa. Such a graph, 
in which the same data as in figure 1 are represented, is illustrated in 
figure 2. The data are on smooth curves which deviate more or less 
from the straight line, 



V(RT~pV) _3 
V c (RT~p c V c ) 2" 



'2 X ' 
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Numerous attempts to find directly a functionfto represent the 
curves in figure 2 have failed — the functions always yielded graphs 
with too much curvature at low densities. 

When the slopes of the curve for C0 2 in figure 2 were plotted against 
density, the resulting graph was remarkably similar to that for the 
function y=ae~ bxn . This suggested the general equation 

V C (R1 C — p c V c ) Jx 

After determining most of the constants / one may write the equa- 
tion in the form 

RT C 

which fulfills the requirement that when X=l, p=p c , bp/bX=0, and 
b 2 p 2 li)X 2 =0, regardless of the value of the constant n; and which 
reduces to 

V(RT e -pV) _S 1 m 

V c (RT c -p c V c )-2 2*' w 

when the quantity RTJp c V c =3. 

When it is desired to calculate pressures corresponding to given 
densities, eq 1 may be expressed in the form 

f.=(^:>--lldft-0+(IS-I)^-™--«] 



+ 



X 3 / BT, 
2 



(Ji-0 



For use in calculations with the equation in either form, eq 1 or 3, 

values of e^ /nH1 ~ xH) dX are given in table 1 for various values of 

1/n and X. For data of moderate accuracy, for example, pressures 
accurate to a few parts in a thousand, graphic methods of interpola- 

i The procedure is as follows: For brevity let F represent temporarily the fraction RTJp t V- Then the 
solution of the general equation for pipe gives 

VlPc=*FX-LX*{F-l)+MX*{F-l) -aX*(F-l)V e-«"dX 

By differentiation 

s ^-F-2LX(F-l)+ZMXKF-l)-2aX(F-l)P e->x n dX+aX2(F-l)e->* n t 

and 

?£t£ e —-2L(F-l)+bMX(F-l)-Sa(F-l) [ l e-*x n dX+4aX(F-l)e->* n -abnX*"(F-l)e-bx\ 

By substitution of the critical conditions p/p.-l, X=l, b(p/p c )l?>X=0, and dHplp e )/dX*=0 in these 
equations, one obtains three simultaneous equations, from which some of the constants can be eliminated. 
The solution is 

. 6-4F/(iM) 

2-M-Fj(F-l) 

Apparently M can have only one value (one-half) if a simple solution is to be obtained. Assumption of 
this value leads to values for the other constants as written in eq.l 
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tion are satisfactory, but for greater accuracy, it may be better to 
calculate directly the values of the integral 2 for the chosen value of n 
and the desired values of X. For values of X up to 1.1 or 1.2, the 
integral is easily calculated by expanding e~^ ,n)xn into series form 
and integrating each significant term of the series. For greater 
values of X this method becomes increasingly laborious, and numerical 
integration of the quantity e ~ Wn)xn becomes the easier method. 
Where data are to be calculated for a large number of values of X 
between 0.8 and 1.2, it is easier to calculate only a few values of 
pressure in this range and to obtain the remainder by interpolation 
between these pressures. 

To maintain the correct relations at the critical point, the quantity 
RTjpcVc should be treated consistently as an exact number in 
calculating the coefficients of the various terms in eq 3, and it is 
advisable in calculating pressures to carry each term of the equation 
to more decimal places than the accuracy of the critical data would 
warrant. 

Table 1. — Values of I e (4/n)(1 ~ xn) dX/or various values of X and 1/n 

JX 



X 


Values of 1/n 





















0.05 


0.1 


0.15 


0.2 


0.25 


1/3 





1.0 


1. 2104 


1. 4428 


1. 6997 


1. 9833 


2. 2965 


2. 8924 


.1 


.9 


1. 0882 


1. 2936 


1. 5175 


1. 7607 


2.0247 


2. 5132 


.2 


.8 


.9661 


1. 1444 


1. 3353 


1. 5382 


1. 7530 


2. 1357 


.3 


.7 


.8439 


.9952 


1. 1531 


1.3158 


1. 4824 


1.7645 


.4 


.6 


.7218 


.8460 


.9709 


1. 0943 


1. 2147 


1. 4065 


.5 


.5 


.5997 


.6969 


.7893 


.8751 


.9541 


1.0710 


.6 


.4 


.4775 


.5479 


.6092 


.6616 


.7063 


.7675 


.7 


.3 


.3554 


.3996 


.4333 


.4591 


.4793 


.5046 


.8 


.2 


.2333 


.2539 


.2668 


.2754 


.2815 


.2878 


.9 


.1 


.1124 


.1166 


.1186 


.1197 


.1204 


.1212 


1.0 























1.1 





-.0696 


-. 0773 


-. 0794 


-.0803 


-. 0809 


-.0814 


1.2 





-.0781 


-. 1080 


-. 1193 


-. 1253 


-. 1286 


-.1317 


1.3 





-.0781 


-. 1124 


-. 1328 


-. 1452 


-. 1526 


-. 1601 


1.4 





-.0781 


-. 1127 


-. 1353 


-.1517 


-. 1627 


-. 1747 


1.5 





-. 0781 


-.1128 


-. 1358 


-. 1532 


-. 1662 


-. 1813 


00 





-. 0781 


-.1128 


-. 1359 


-. 1534 


-. 1673 


-. 1855 



and when X>1, 
When n=0 



The constant n is the only quantity in eq 1 which is not explicitly 
determined by the critical data. The most convenient method for 

* It may be of academic interest to note that for n= <*>, when X<1, 

\\w>M-z n )dX=l-X 

fJg(</«)(i-x")dX=0. 

As an aid in checking calculations, it may be of interest to note that 

J^V»"<fy=(l/70r(l/7i). 

See Bierens de Haan, Nouvelles Tables d'Integrales D6finies, 1867, table 81, formula 8. (Printed by J. C. 
Drabbe.) Formulas for calculating (l/n)T(l/n) are given in Jahnke-Emde Funkionentafeln 3d (revised) 
edition, p. 10, published by G. B. Teubner, Leipzig and Berlin (1938). Then 

J^ OT eWn)(i-x-)d^=( n /4) , 'V / "J^ 00 e-»"dy=(w/4)i/»«eV*(i/n)r(l/n). 
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determining n is to substitute the critical data in the equation and to 
calculate the quantity 

/ V(RT c -yV) 3,l y \//3 1 \ 

\V c (RT-p c V c ) 2 i "2 A //V2 l-{ Vc V c )IRTj 

which is equal to 

f e WnH1 ~ xn >dX. 
Jx 

The corresponding value of n is obtained by interpolation in table 1 . 
For this purpose, experimental data at slightly less than one-half 
the critical density are most desirable, since calculations from eq 1 at 
various densities with the same critical data, but with two slightly 
different values of n, yield two sets of calculated pressures which 
differ by a maximum amount at a little less than one-half the critical 
density and which become equal as X approaches or 1. 

Of the three critical data, the critical volume is usually know r n with 
the least accuracy. Fortunately, however, there is a compensating 
action between the various terms of eq 1, such that a change in the 
critical volume alone produces very little change in the pressures 
calculated at given densities. 

A limit which the constant n can neither equal nor exceed, if eq 1 
is to represent properly the critical state, is afforded by consideration 
of the third and higher derivatives of pressure with respect to density. 

The value at the critical state of the third and higher derivatives is 
probably beyond experimental determination. It is reasonable to 
suppose, however, that for densities less than the critical density, the 
pressure does not exceed the critical, and that for densities greater 
than the critical density no pressure less than the critical exists. This 
supposition leads to the conclusion that if all derivatives of pressure 
with respect to density up to and including any derivative of odd order 
are zero, then the derivative of next higher (even) order must also be 
zero. 

One finds by differentiation of eq 1 for the condition X=l, that 
ZPp/dX 3 is zero when 



K#-)/(iH) 



and that within the useful range of n, d 4 £>/dX 4 is zero only when n is 
approximately 6.7. Hence one may say that in general n can neither 
equal nor exceed the ratio 



(#:-)/(iH> 



The fact that the values of the constant n determined for various 
substances usually are less than 



(#-0/(iH) 



by considerable amounts does not preclude the possibility of d 3 p/dX 3 
(and d 4 2?/dX 4 ) actually being zero at the critical point. The equation 
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is only empirical and therefore not exact. In the neighborhood of the 
critical point only very minute differences exist between the pressures 
calculated from the values of n which were used and from those which 
equal 



(^i-O/dft-0- 



For example, the reduced pressure of normal heptane calculated 
for X=0.9 and n = 3.66 is 0.9997, whereas the use of n=5 or 
approximately 



(tifrWte) 



gives 0.99997 for the reduced pressure — a difference of about 3 parts 
in 10,000. The experimental data in this case indicate an inter- 
mediate value. 

The method described for determining n from experimental data 
was used for H 2 0, (RT e lp e V e =4.Z35, n=3.12); NH Z (RTJp c V c = 
4.075, n=3.20) and for C0 2 (RT c lp c V c =3M9, n=5A8). A value 
of n for H 2 (i?T c /^ c T 7 c =3.203, n=16.6) was chosen more or less arbi- 
trarily as being equal to 



(*$H)/' 



since the experimental data were best represented by the use of as 
large a value of n as possible. Although this procedure is mathe- 
matically incorrect, the errors introduced are practically negligible. 
As indicated in table 3, the calculated reduced pressure at 3l==1.I is 
0.9999, that is low by 1 or more parts in 10,000. 

In figure 3 the value of RT c /p c V c is plotted as ordinate and that of 
\\n as abscissa. The dotted curve indicates the values of \jn for 
which b 2 p/dX z =0 at the critical state, whereas the values for the 
substances mentioned are indicated by the labeled crosses. The 
continuous, straight line in figure 3 through the point for H 2 and that 
for C0 2 passes between the points representing water and ammonia. 
The values of n given in table 2 for other substances were obtained 
from this line and were found to represent the data within the experi- 
mental accuracy. 

IV. REPRESENTATION OF THE EXPERIMENTAL DATA 

Although table 2 contains the critical temperatures, pressures, and 
volumes for the substances discussed in this paper, the main purpose 
of this table is to present the empirical constants actually used in 
representing the critical isotherms. For this reason some of the data 
are given to more places than is warranted by the experimental 
accuracy. As noted at the bottom of the table, data differing slightly 
from those chosen by the authors have been used in some instances. 
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Figure 3. — The value of ljn as a function of RT e lp c V c , 
Table 2. — Constants used in equations for the critical isotherms 



Substance 


T e * 


P e 


Vc* 


PcVc 


P e V e 


1/71 


Refer- 
ence 


Hydrogen 


°K 

33.18 

154. 27 
125. 96 
191.1 
305. 43 
369. 96 
f 304. 197 

'~~364.~56~~ 

425. 16 
466.6 

540. 17 
j 406. 1 


atm 
f (13. 03) 
\ 12. 75 
49. 713 
33. 490 
45.8 
48.20 
42. 004 


0. 002909A 


3.203 
3.273 
3.422 
3.422 
3.448 
3.513 
3.592 
3.649 


0. 3122 
.3055 
.2922 
.2922 
.2901 
.2847 
.2784 
.2740 


0. 0601 
.080/ 
.120 
.120 
.128 
.145 
.167 
.1823 


1 






Oxygen 


. 003324A 
. 004027A 
.0993 
.148 
(• 201) 
. 004216A 


1 


Nitrogen 


1 


Methane.. . 


2 


Ethane 


3 


Propane 


4 


COj_ 


5 


72.885 
45.4 
37. 468 
36.30 
27.00 






6 


Propylene 


.180 
(.2544) 
.279 
.416 


3.660 
3.660 
3.774 
3.946 


.2732 
.2732 
.2650 
.2534 


.184 
.184 
.218 
.264 


7 


n-Butane 


8 


Ethyl ether 


9 


n-Heptane 


10 




11 


Ammonia 


112.1 










12 




fc::::~:: 


.0730 


4.075 


.2454 


.312 


13 




f 647. 27 


219.4 


14 




.0558 


4.335 


.2368 


.32 


15 











» The temperatures for H2, O2, and N2 are as reported in the reference cited. In other instances, °K= 
273.16-f°C. 

b Liters per mole, unless marked "A", in which cases Amagat units were used. The numbers in paren- 
theses are different from those given in the references cited. The reasons for these changes are explained in 
the text. 
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A knowledge of the weight of the normal liter is involved in the use 
of Amagat units, but this quantity is so much more accurately known 
than the other data involved that it is considered unnecessary to 
quote the source of the data used for the normal liter. 

1. HYDROGEN 

Since a discrepancy of about 2 percent exists between the experi- 
mental data for the superheated vapor of hydrogen and those for the 
saturated vapor, 3 two sets of constants were used in eq 3 for calculating 
values of p/p e - One set with p c = 13.03 atm (see table 2) was used in 
calculating values used for comparing with the superheat data in 
figure 4, and a second set was used with p c — 12.75 atm for com- 
parison with tables which may be based upon this value of the critical 
pressure. The values calculated from these two sets of constants are 
given later in columns 2 and 3, respectively, of table 3. Although 
there is a considerable difference at low densities between the reduced 
pressures calculated from these two sets of constants, the absolute 
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Figure 4. — Comparison of pressures calculated from equation 3 with pressures 

observed for hydrogen, 

+ Interpolated from data by G. P. Nijhoff and W. H. Keesom in Comm. Phys. Lab. Univ. Leiden No. 

188e, 17 (1924-28). 
O Averages of data by H. K. Onnes and P. M. Penning in Comm. Phys. Lab. Univ. Leiden No. 165b, 

15 (1920-26). 
X C A. Crommelin and^J. C. Swallow,'Coinm. Phys. |Lab. ;Univ.^ Leiden No. 172a 16 (1924-26). 

values of pressure calculated for X=0.1*or*.less differ]jby~onlyj r[or^2 
parts in 1,000. 

Figure 4 presents a comparison of experimental and calculated 
pressures for hydrogen. The agreement is as good as can be expected 
in view of the discrepancy already mentioned. The two plus marks 
at X= 0.058 represent the largest and the smallest pressures which 
were obtained by different assumptions in interpolating between the 
experimental data at other temperatures. The trend of the crosses 
indicates that this isotherm is above the critical temperature. This 
fact, combined with the 2-percent discrepancy already mentioned, 
suggests a possible error in temperature measurement. 



2. OXYGEN AND NITROGEN 

The experimental data for both oxygen and nitrogen are at irregular 
intervals and therefore not subject to easy interpolation to the critical 
temperature. For the purpose of comparing these data, use is made 
of an equation developed by Leib, 4 which, for diatomic gases, may be 
written in the form 

pV_ K(T c /Typy c /v 

RJ 1 e K(T e ITl)V>V e IV_l 



» H. K. Onnes, C. A. Crommelin, and P. G. Cath, Comm. Phys. Lab. Univ. Leiden No. 151c, 14 (1914-22). 
«E. F. Leib, Trans. Am. Soc.JVIech. Engrs.^63, 167.(1941). 
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and may be expanded into the form 

^F/i?r=i-(i^/2)(r c /r) 5 / 2 ^F+(K(r^T) 5 /HyF)7i2- .... 

At low densities this leads to 

mm m ivm a -pv/rt) + (K/2y(T c /Ty< 2 (v c /v)/z=K/2, 

where K/2 is a single arbitrary constant which may be determined 
from the experimental data. Equation 1 may be expanded to a some- 
what similar series, in which 1.157 corresponds to K/2 and 0.454 to 
the third virial coefficient (see table 4). According to Leib, the third 
virial coefficient would then be one-third of (K/2) 2 , or 0.446. With 
the choice of coordinates in figures 5 and 6, the left-hand portion of 
the continuous curve representing eq 1 must be horizontal. The 
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Figure 5. — Comparison of data calculated from equation 1 with data observed for 

oxygen. 



-154.27°K. Calculated from equation 1. 

-154.27°K. From E. F. Leib's equation, Trans. Am. Soc. Mech. Engrs., 63, 157 (1941). 
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Onnes and H. A. Kuypers, Comm. Phys. Lab. Univ. Leiden No. 169a, 16 



X 

? 

dotted curve just above the continuous curve represents Leib's equa- 
tion. Leib's theory indicates that with the introduction of (T c /T) 5 / 2 , 
the data for the different isotherms should lie on a single curve, which 
is approximately the case at the lower densities. Figure 5 indicates 
that eq 1 represents the data for oxygen within the precision of the 
experimental data. The data for nitrogen are more precise, and it 
would appear in figure 6 that at the higher densities, the critical iso- 
therm calculated from eq 1 is slightly too low to correspond to the 
uniform spacing per degree of the different isotherms. This must be 
interpreted as a discrepancy between the reported critical data and 
the other isothermal data, rather than a failure of eq 1, since the form 
of the calculated isotherm appears to be correct and, at the critical 
state, the equation is determined exactly by the critical data substi- 
tuted in it. 

The dotted curve at the bottom of figures 5 and 6 provides means 
for estimating the errors in pressure equivalent to deviations in the 
ordinate. 
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All observations made on the oxygen, using a piezometer with a 
reported normal volume of 1820.6, 5 are designated in figure 5 with a 
prime mark. These points deviate systematically from smooth curves 
through the other data. The deviation corresponds in size to an 
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Figure 6. — Comparison of data calculated from equation 1 with data observed for 

nitrogen. 



125.96° K. Calculated from equation 1. 

125.96° K. From E. F. Leib's equation, Trans. Am. Soc. Mech. Engrs., 63, 157 (1941). 
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error of 1 gram in the weight of mercury in calibrating the small end 
of that particular piezometer. 

3. METHANE 

There are no data for methane at temperatures near the critical, 
but the calculated pressures are in agreement with those extrapolated 
from the observed data. 6 

4. ETHANE 

A comparison of the pressures calculated from eq 3 with the observed 
data for ethane is given in figure 7. The critical constants were those 
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Figure 7. — Comparison of pressures calculated from equation 3 with pressures 

observed for ethane. 

X J. A. Beattie, C. Hadlock, and N. Poffenberger, J. Chem. Phys. 3, 93; 1935. 

O J. A. Beattie, Gouq-Jen Su, and G. L. Simard, J. Am. Chem. Soc. 61, 924 (1939). 

obtained from the data near the critical density, represented by 
circles in the figure. The observed pressures at lower densities repre- 

* H. K. Onnes and H. A. Kuypers, Comm. Phys. Lab. Univ. Leiden No. 169a, 16, 5 (1924-26). 
e H. M. Kvalnes and V. L. Gaddy, J. Am. Chem. Soc., 53, 394 (1931). 
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sented by crosses are lower than the calculated values. This differ- 
ence appears to be due to a difference between the two sets of experi- 
mental data rather than a failure of the equation. Possibly there is a 
difference in the purity of the two samples. 

5, PROPANE 

The^critical volume of propane 0.201 liter per mole was used in 
place of the value 0.195 chosen by Beattie, Poffenberger, and Had- 
lock. The value used here corresponds to the mean of the densities 
of the vapor and of the liquid at 96.70° C (about 0.1° below the 
critical temperature). The pressures calculated from eq 3 are com- 
pared with the values observed for propane in figure 8. Values 
calculated from the equation agree with the data represented by 
circles, and agree as well as can be expected with the data for the 
other sample (represented by crosses), since there is a consistent 
difference of about 1 part in 1,000 between corresponding pressures 
for the two samples. 

6. CARBON DIOXIDE 

Two previous papers in which the critical temperature was chosen 
first 7 near 31.1°C and later 8 at 30.96°C have appeared to be rather 
contradictory. This apparent contradiction seems to be cleared up 
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-Comparison of pressures calculated from equation 3 with pressures 
observed for propane. 

X J. A. Beattie, W. C.'Kay, and J. Kaminsky, J. Am. Chem. Soc/59,^1589 (1937). 
O^J. A. Beattie, N. Poffenberger, and C. Hadlock, J. Chem. Phys. 3,J)6 (1935). 

by the theoretical work of Mayer and Harrison 9 , in which they 
conclude that a temperature exists for which, at the critical density, 
P—Pc, and the first and second derivatives of pressure with respect 
to volume are zero; but that the sharp dividing surface between the 
liquid and the gas phases disappears at a slightly lower temperature. 
For the present purpose, 31.037°C was used as the critical temperature, 
that being the value reported by Michels. The critical density is 
chosen from Michels' data, but the pressure is that observed at the 
National Bureau of Standards. 

The agreement between the two sets of data (see fig. 9) is remark- 
able in view of the fact that volumes were only roughly measured in 
the latter investigation. The calculated pressures agree very well 
with the observed data, as is shown in figure 9, the maximum deviation 
only slightly exceeding 3 parts in 10,000. There is a possibility that 
these small deviations might be further reduced if the data had been 
at a temperature about 0.05°C higher. 

* C. H. Meyers and M. S. Van Dusen, BS J. Research 10, 381 (1933) EP 538. 
» H. T. Kennedy and C. H. Meyers, Refrig. Eng. 15, 125 (1928). 

• J. E. Mayer and S. F. Harrison, J. Chem. Phys. 6, 87 (1938). 
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The pressures calculated from eq 3 in the range X=\.l to 1.3 
are apparently too large. As the density is increased beyond X=\.?>, 
the calculated pressures soon become too small and rapidly diverge 
from the observed values. For example, when X=\.2 the pressure 
observed at this Bureau is 72.97 atm, the calculated 73.13, or 22 
parts in 10,000 higher. For X— 1.3, the observed pressure is 73.61 
atm, whereas the calculated value is 74.04, or 50 parts in 10,000 
higher. Michels reports 81.261 atm at a density corresponding to 
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Figure 9. — Comparison of pressures calculated from equation 3 with pressures 
observed for carbon dioxide. 

X C. H. Meyers and M. S. Van Dusen, BS J. Research 10, 381 (1933) RP 538. 
0+ A. Michels and C. Michels, Proc. Roy. Soc. (London) [A] 153, 201 (1935). 

X— 1.485. Here, however, the calculated value, 80.08 atm, is 161 
parts in 10,000 lower. In the region where X is greater than 1.1, 
the equation deviates in a similar manner in the case of other sub- 
stances. 

7. PROPYLENE 

The data reported for propylene consist of values of pV/RT tabu- 
lated for even reduced pressures. In figure 10 these values are com- 
pared with those calculated from eq 3. 

Since the comparison is on the basis of pV/RT, a decrease of about 
1 percent in the assumed critical volume would be sufficient to cause 
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Figxjee 10. — Comparison of data calculated from equation 3 with data observed for 

propylene. 
W. E. Vaughan and N. R. Graves, Ind. Eng. Chem. 32, 1252 (1940). 

practical agreement between calculated and observed values. An 
accuracy of only 1 percent is estimated by the authors for the data. 



8. NORMAL BUTANE 



The use in this paper of the critical volume 0.2544 liter per mole, 
together with the other constants in table 2, leads to a value 0.2732 
for p c V c /RT c . The pressures calculated with this value of the critical 
volume, together with n from figure 3, are in slightly better agreement 
with the observed pressures than those calculated with the aid of the 
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reported value, 0.258 liter per mole. The excessive number of decimal 
places in the value 0.2544 was used for convenience to make the 
constants in eq 3 the same as for propylene. 

As shown in figure 11, the pressures observed near the critical state 
and the values calculated from eq 3 are in excellent agreement. The 
difference between the calculated and observed pressures at lower 
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Figure 11. — Comparison of pressures calculated from equation 3 with pressures 
observed for normal butane. 

X J- A. Beattie, G. L. Simard, and Gouq-Jen Su, J. Am. Chem. Soc. 61, 24 (1939). 
O J- A. Beattie, G. L. Simard, and Gouq-Jen Su, J. Am. Chem. Soc. 61, 26 (1939). 

densities amounts to about 1 part in 1,000, but this difference is only 
about one-fifth the difference between the observed pressures and 
values calculated from the equation given by the authors cited in 
figure 11. 

9. ETHYL ETHER 

As illustrated in figure 12, the difference between the pressures 
observed for ether and those calculated from eq 3 is less than 2 parts 
in 1,000. 

10. NORMAL HEPTANE 

The comparison of data on heptane in figure 13 suggests a systematic 
difference between the two sets of observed data. However, this ap- 
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Figure 12. — Comparison of pressures calculated from equation 3 with pressures 

observed for ethyl ether. 
J. A. Beattie, J. Am. Chem. Soc. 46, 342 (1924). 

parent difference may be in part caused by errors in extrapolation 
from 275° to 267° C, for the data represented by crosses. 

11. AMMONIA 

The comparison of data for ammonia in figure 14 shows the agree- 
ment for low densities (see upper portion of figure) to be excellent not 
only between the calculated and the observed pressures but also be- 
tween the two sets of observed data. The deviations from eq 3 
rarely exceed 5 parts in 10,000. A study of a graph similar to figure 2, 
in which other isotherms than the critical are plotted for ammonia, 
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gives the impression that this good agreement between the calculated 
and observed values is somewhat fortuitous. The isotherms at higher 
temperatures bend upward at low densities too sharply to be rep- 
resented by eq 3, whereas the isotherms at lower temperatures are 
similar to those for water in being convex upward. The calculated 
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Figure 13. — Comparison of pressures calculated from equation 3 with pressures 
observed for normal heptane. 

X L. B. Smith, J. A. Beattie, and W. C. Kay, J. Am. Chem. Soc. 59, 1587 (1937). 
O J- A. Beattie and W. C. Kay, J. Am. Chem. Soc. 59, 1586 (1937). 

pressures at higher densities (see lower portion of fig. 14) show larger 
differences from the observed values; but below the critical density, the 
differences are probably due to experimental error and uncertainty 
of extrapolation — in this region there are no corresponding deviations 
of this size for other substances. The differences at densities greater 
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Figure 14. — Comparison of pressures calculated from equation 3 with pressures 

observed for ammonia. 

X C. H. Meyers and R. S. Jessup, Refrig. Eng. 11, 345 (1925). 

O J. A. Beattie and C. K. Lawrence, J. Am. Chem. Soc. 52,6 (1930). 

+ F. G. Keyes, J. Am. Chem. Soc. 53, 9G5 (1931). 

than the critical are mainly due to inadequacy of eq 3, as mentioned 
in the discussion of carbon dioxide. 



12. WATER 

Figure 2 shows that eq 3 does not represent the data for water as 
satisfactorily as it does the data for other substances. The differences 
of over % percent between the observed and calculated pressures, 
illustrated in figure 15, appear to be greater than can be attributed to 
experimental error. Water differs from ammonia in that its isotherms 
when plotted on a graph similar to figure 2 are all convex upward at 
low densities. 
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Table 3 contains values of reduced pressure calculated from eq 3 for 
13 substances at various reduced densities, using the constants given 
in table 2. In some instances the value of the integral in the equation 
may not have been calculated with sufficient accuracy to determine the 
corresponding value of the reduced pressure to better than two or 
three units in the last decimal place. The values of reduced pressure 
for densities up to and including X=\.\ are believed to agree with the 
observed data within reasonable A limits of experimental accuracy 
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Figure 15. — Comparison of pressures calculated from equation 8 with pressures 

observed for water. 

F. G. Keyes, L. B. Smith, and H. T. Gerry, Proc. Am. Acad. 70,319 (1935). 

except in the case of water, where some of the calculated values appear 
to be in error by about % percent. 

Table 3. — Reduced densities and corresponding reduced pressures calculated from 
equation (8) for various substances at the critical temperature 





Reduced pressures 


X 


•Hj 


fcHa 


Oaand 
Ni 


CH* 


C 2 H(j 


C 3 H 8 


COj 


C 3 H« 
and 


Ether 


C 7 Hu 


NH 3 


H 2 
















C4H10 










0.1 


0. 2872 


0.2927 


0. 3041 


0. 3061 


0.3110 


0. 3169 


0. 3211 


0. 3219 


0. 3300 


0. 3419 


0. 3500 


0. 3689 


.2 


.5132 


.5215 


.5384 


.5413 


.5484 


.5566 


.5624 


.5636 


.5745 


.5898 


.5991 


.6249 


.3 


.6854 


.6943 


.7122 


.7151 


.7223 


.7305 


.7361 


.7373 


.7472 


.7606 


.7670 


.7920 


.4 


.8109 


.8190 


.8347 


.8372 


.8432 


.8497 


.8540 


.8550 


.8620 


.8711 


.8735 


.8935 


.5 


.8992 


.9038 


.9154 


.9170 


.9212 


.9255 


.9281 


.9288 


.9325 


.9368 


.9364 


.9496 


.6 


.9522 


.9563 


.9634 


.9643 


.9664 


.9685 


.9699 


.9702 


.9716 


.9737 


.9708 


.9779 


.7 


.9826 


.9846 


.9878 


.9881 


.9890 


.9896 


.9900 


.9901 


.9901 


.9900 


.9884 


.9914 


.8 


.9960 


.9967 


.9975 


.9975 


.9977 


.9977 


.9977 


.9977 


.9975 


.9972 


.9966 


.9977 


.9 


.9997 


.9998 


.9998 


.9998 


.9998 


.9998 


.9998 


.9998 


.9997 


.9996 


.9995 


.9997 


1.0 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1.1 


•. 9999 


1. 0000 


1.0000 


1.0002 


1.0002 


1.0003 


1. 0003 


1. 0003 


1.0005 


1.0004 


1. 0005 


1. 0005 


1.2 


1.0012 


1.0016 


1.0015 


1. 0022 


1. 0025 


1.0030 


1. 0033 


1.0035 


1. 0041 


1.0051 


1.0062 


1.0060 


1.3 


1. 0141 


1. 0130 


1.0127 


1. 0132 


1.0138 


1.0150 


1.0158 


1.0161 


1.0181 


1. 0210 


1. 0241 


1.0241 


RT C 

PeV, 


3.203 


3.273 


3.422 


3.448 


3.513 


3.592 


3.649 


3.660 


3.774 


3.946 


4.075 


4.335 



• Calculated with p c = 13.03 atm, l/n=0.060 (see table 2). 

* Calculated with p e = 12.75 atm, l/rc=0.080 (see table 2). 

« Actually calculated as less than unity because of choice of constants (see text). 

As already mentioned, two sets of values are given for hydrogen, 
since the choice of a value for the critical pressure of that substance 
is to some extent a matter of personal judgment. 

The value of RT c JpcV e is also given at the base of the table as an 
aid in comparing with data for other substances not included in the 
table. For example, the ratio RTJp c Y c for argon is very nearly the 
same as that for oxygen and nitrogen, consequently the reduced 
pressures given for those two substances should apply also to argon. 
It may be convenient also to represent the reduced pressures for other 
substances by interpolation between the values given in table 3. 
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V. VIRIAL COEFFICIENTS FOR THE ISOTHERM AT THE 
CRITICAL TEMPERATURE 

Equation 1 may be transformed to 



2V RTJ 

Or approximately for small values of X 

■^^'-SXK'-S')- 1 )!''™- 1 ""] 



(4) 



~RT~ 



-X 2 



(5) 



The range over which this approximation fulfills a given requirement 
depends upon the value of n. For example, pV/RT calculated from 
eq 5 does not differ more than 1 in 10,000 from the same quantity 
calculated from eq 4 for values of X up to 0.5 or 0.6 for H 2 ; up to 0.3 
or 0.4 for C0 2 ; and up to 0.1 for NH 3 . 
One may rewrite eq 5 briefly as 



pV/RT=l+BX+CX 2 . 



(6) 



Table 4 gives values of B and G for various substances as calculated 
from eq 5. The virial coefficients, as sometimes defined from the 
equation 

can be derived easily from these coefficients. 

Table 4. — Virial coefficients for critical isotherms 





Coefficients from 






Substance 


Equation 5 


Berthelot 






-B 


C 


-B 


-B 


C 


Hydrogen _ 


1.072 
1.157 
1.157 
1.164 
1.195 
1.230 
1.256 
1. 261 
1.261 
1.318 
1.413 
1.492 
1.548 


0.384 
.454 
.454 
.460 
.488 
.521 
.548 
.552 
.552 
.613 
.729 
.837 
.902 


1.126 
1.203 
1.203 
1.212 
1.235 
1.263 
1.283 
1.287 
1.287 
1.333 
1.387 
1.433 
1.524 


» 1.05 

hi. 23 
« 1.23 




Oxygen _ 


b 0. 595 


Nitrogen.. _ __ 


c.586 


Methane _ 




Ethane _ 






Propane _ 






Carbon dioxide__ _. 


d 1.261 


<*.580 


Propylene 




n-Butane 






Ethyl ether 






w-Heptane _ 






Ammonia 






Water* 













* G. P. Nijhoff and W. H. Keesom, Comm. Phys. Lab. Univ. Leiden No. 188e 17 (1924-28). 
b H. A. Kuypers, Comm. Phys. Lab. Univ. Leiden No. 169b, 16 (1924-26). 

o A. Th. van Urk, Comm. Phys. Lab. Univ. Leiden No. 169e, 16 (1924-26). 
a Proc. Roy. Soc. (London) [A] 153, 201 (1935). 

• Allowing for the effect of C, the value of B estimated from the 2 lowest observed pressures is about 1.46. 
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The fourth column of table 4 contains values oi—B calculated from 
Berthelot's equation 



p v = R i i +MM 1 ~ 6 ¥)] 



which leads to a value — 45RT c /12Sp e V c for B at the critical tempera- 
ture. Comparison values presented by other authors are given in some 
cases. 

In general, the agreement of pressures calculated from eq 3 has been 
shown to be very good, hence the values of B in table 4 obtained from 
eq 5, with the exception of that for water, are considered to be within 

2 percent or less of the correct values. 

The accuracy of the data on hydrogen is not a sufficient basis for a 
choice between the coefficient calculated from eq 5 and that given by 
the Leiden Laboratory. 

The data on oxygen and nitrogen confirm the coefficients presented 
in table 4 from eq 5. The value 1.23 given by the Leiden Laboratory 
for B for oxygen and nitrogen, in other words, the intercept for the 
critical isotherm at zero density in figures 5 and 6 is out of line with the 
experimental data for both substances. 

The discrepancy of about % percent between the value of B for C0 2 
calculated from eq 5 and that given by Michels, is insignificant for the 
purpose of calculating pressures at low densities. For example, for 
1 atm eq 5 gives pV/RT=0. 99526, whereas Michel's equation gives 
0.99524; even at 19 atm the difference in pV/RT amounts to only 3 
parts in 10,000. 

It will be seen from table 4 that the values of B from Berthelot's 
equation are about 5 percent too high at the top of the table, but for 
larger values of RT c /p c V c (see table 2), the difference between the 
two sets of values decreases and changes sign between ethyl ether and 
n-heptane. 

The deviations shown in figure 15 indicate that pressures calcu- 
lated from eq 3 are either in agreement with or slightly greater than 
the observed data for ammonia at the lowest densities. Consequently, 
the value of B from eq 5 in table 4 is not too large, and the value of B 
from Berthelot's equation is about 4 percent too small. 

The fact that eq 5 and Berthelot's equation lead to values of B for 
water respectively 6 and 4 percent too large, indicates that water is an 
exceptional substance, in fact, the only one encountered which eq 5 
obviously does not represent within the limits of experimental accuracy. 

From Leib's equation already mentioned in the discussion of oxygen 
and nitrogen, the third virial coefficient C is equal to 5 2 /3, a ratio 
which is in agreement with eq 5 only when the value of RT c /p c V e is 

3 for the substance considered. For substances having larger values 
of RT c /p e V c , smaller values of C are calculated from Leib's equation 
than from eq 5. The larger is the value of RT c /p c V Cy the greater is 
the difference between the two values of C. For ammonia the 
difference between the two values of C is about 11 percent. 
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VI. FUGACITY ALONG THE ISOTHERM AT THE CRITICAL 

TEMPERATURE 

For the purpose of calculating fugacity along the critical isotherm, 
an expansion of the integral appearing in eq 3 into series form is con- 
venient. Then solving for 9 (p/p c ) /dX by differentiating that equation, 
substituting this derivative in the equation 

and integrating, one obtains 

fi 3p c FA /w f 4/n (-l)"(4/n)» 1 

V 1 BT.J 6 i 1 n+l + ■ ' ■ f (mft+l)(m!) + ".'if 



"•"2V ST./ J2 A (»+l)(n+2) + 

(1— ) w (4/7i)" 1 (mu+3)Z""' +2 

+ (m»+l)(m»+2)(m!) + * ' ' ' 



(7) 



For a hypothetical substance with a value of 3 for RT c /p c V c , eq 7 
becomes 

>^>°#-<'-»-f> 

and for the critical state of such a substance, \n(f c /p c ) = (hi 3) — 1.5* 
or f c /Pc= 0.6694. For carbon dioxide with the constants in table 2, 
c/Pc=0.Q593. This is in perfect agreement with the value inter- 
polated between the data given by Deming and Deming. 10 

For carbon dioxide under a pressure of 1 atm, the equation presented 
in this paper leads to pV/RT=0.9952Q and to//p = 0.99528, whereas 
the corresponding value given by Deming and Deming for f/p is 
0.99504. The fact that Michels obtains pV/RT=0.99524: and that 
i/p should approximate pV/RT very closely, confirms the value 0.99528 
or f/p for carbon dioxide at 1 atm and the critical temperature. 

It is to be noted that the value of f/p at the critical state is very 
slightly over 2/3 when RT c /p c V c =3 and decreases only slightly with 
increase in RT c /p c V c , so that for any substance at the critical state, 
f/p deviates from 2/3 by only a few percent. 

VII. SUMMARY 

The equation 

V(RT c -pV) _3 1/3 1 \ r e ^ n) a-xn )dx 

V c (RT-p c V c ) 2 2 Ai A2 l-PcVJRTjJ x e a *> 

where n is an empirical constant and X—V c /V, the reduced density, 
has been shown to fulfill the requirements at the critical state and to 

i° W. E. Deming and L. S. Deming, Phys. Rev. 56, 108 (1939). 
469958—42 5 



176 Journal of Research of the National Bureau of Standards 

represent the critical isotherms for 12 out of 13 substances in the range 
of densities up to 1.1 times the critical density, within the reasonable 
limits of experimental error. The deviations for the one exception, 
water, amount to % percent, which exceeds the tolerance set by the 
experimental accuracy. For the very accurately measured critical 
isotherm of carbon dioxide, the differences between calculated and 
observed pressures do not exceed 3 parts in 10,000. 

Of the five constants in the equation, one is limited to a single 
value for a simple solution, three are directly determined by the critical 
data, which when substituted in the equation make it automatically 
fulfill the requirement of the critical state, that the first and second 
derivatives of pressure with respect to density shall be zero. The 
reciprocal of the fifth constant, n, has been shown to be very nearly a 
linear function of RTJp c V c . In fact, when RTJp c V c does not exceed 
3.8 or 3.9, this linear relation established by the values determined for 
hydrogen and for carbon dioxide furnishes a correlation from which n 
may be determined for other substances, and consequently, the critical 
isotherm may be predicted with great accuracy for densities up to 1.1 
times the critical density. The value of n for substances with larger 
values of RT c /p e V c will probably need to be determined individually, 
although even then pressures may be predicted within a few percent, 
from the critical data and the linear relation mentioned for the fifth 
constant. 

For low pressures the equation may be reduced to one expressing 
pV/RT as a quadratic function of the density, from which pressures 
may be calculated with an accuracy exceeding that of most measure- 
ments. 

Although the equation is applied only to the critical temperature, 
the wide range of densities covered with such accuracy makes it valu- 
able for establishing the critical isotherm as a base line in the tabula- 
tion of smoothed values of the properties of substances. 

Values of fugacity along the critical isotherm can be calculated from 
the equation. The fugacity at the critical state is very nearly two- 
thirds the critical pressure for all substances, but varies slightly with 
the value of RT c /p e V e . 
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